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RHIC data on high-p^ hadron production display strong suppression in a wide rapidity 
region, indicating strong induced energy loss for both transversally and longitudinally 
traveling partons. We investigate the interplay of energy loss and rapidity dependence in 
a perturbative QCD improved parton model, and estimate the opacity of the produced hot 
p • matter in AuAu collisions at energies ^/s = 200 AGeV and 63 AGeV at different rapidity 
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values. Direction-dependent suppression offers the possibility to study the geometry of 
D . the hot matter. 



1. Introduction 

>< 

^ I In heavy-ion collisions at RHIC energies we expect the formation of hot and dense mat- 
ter in the central space-time region of the collision. One can measure the density profile of 
the color charges in this deconfined region using induced energy loss of partons traversing 
the dense region in different directions. Data from central rapidity show the transversally 
propagating partons and their large energy loss (see AuAu data from PHENIX [ 1 and 
STAR [|2j). At large rapidities strong suppression has been seen by BRAHMS [|21|^ in 
pion and hadron production. This finding offers the possibility to analyze energy loss for 
different momentum directions. 

We apply a perturbative QCD improved parton model to describe hadron production 
at high pt [EllEj- The Glauber-Gribov model provides the proper framework to describe 
hadron production as a superposition of proton-proton results to obtain baseline hadron 
spectra in heavy-ion collisions. The parton energy loss is calculated by the GLV-method [ 
|7j, and the opacity of the dense matter is extracted in different rapidity regions, i.e. for 
different parton momentum directions. We include shadowing and multiscattering, tested 
in a wide rapidity range in dAu collisions [ |H1 El • 
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2. Theoretical Background of the Model 

Pion production is calculated in AuAu collisions at RHIC energies using a perturbative 
QCD improved parton model in leading order (LO) [El El- Schematically 

E^—^r^Y.TAA'{r,h)0fa/piXa,'kTa,Q^)®fb/A{Xb,'kTb,Q^)®da0D^{z„Q^) . (1) 
P abc 

The collision geometry and the superposition of the nucleon-nucleon collisions is in- 
cluded by the Glauber-Gribov model. In the multidimensional convolution integral 
TAyi'(r, b) = t^(r)tyi/(|b — r|) is the nuclear overlap function for the colliding nuclei. The 
nuclear thickness function for nucleus A is introduced as tA^b) = J dz pA^b, z), normalized 
as / d'^btA{b) = A. We use a Woods- Saxon distribution for Au. 

Nuclear PDFs can be constructed from proton and neutron PDFs [E]: 



fa/A{x, Q^) = Sa/A{x) -fa/p{x, Q^) + - ^) fa/n{x, Q 



(2) 



The function Sa/A{x) describes 'conventional' nuclear shadowing. Here a fe-independent 
Sa/A{x) was taken from the HIJING parametrization [llOj. 

The last term in eq. (Q) is the fragmentation function (FF), D'^{zc,Q'^). This term 
gives the probability for parton c to fragment into a pion with momentum fraction Zc at 
scale Q = kpjn. We apply the KKP parametrization [E]- The factorization scale [Q) is 
connected to the momentum of the intermediate jet, Q = k, ■ px / z^, which scales together 
with Q in the FF. We applied k = k = 2/3 at RHIC energies. 

The density of the color particles can be measured by the induced energy loss of high- 
energy quarks and gluons traveling through the hot dense deconfined matter. One can 
determine the non-abelian radiative energy loss, AE{E, L). This quantity depends on the 
color charge density through the opacity, n = L/X, where L is the length of the traversed 
matter and A is the mean free path. In "thin plasma" approximation the energy loss is 
given by the following expression [[7]: 

Here Cr is the color Casimir of the jet and yU^/ A ~ a^Ppart is a transport coefficient of the 
medium, which is proportional to the parton density, Ppart- The color Debye screening 
scale is denoted by p. N{E) is an energy dependent numerical factor with an asymptotic 
value 4 at high jet energies. 

Jet energy loss influences the flnal hadron spectra. We can introduce this effect via 
modifying the momentum fraction of the outgoing parton before the fragmentation. Con- 
sidering an average energy loss, AE, in a static plasma, the argument of the fragmentation 
functions will be modifled as 



Dn/c{z,,Q^) _^ 4 D^/,{Z*,Q'') 
-nzl Zc T^zl 



(4) 



Here z* = z^j (1 — AE jp^) is the modifled momentum fraction. 
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3. Nuclear modification factor, -R^^, at large rapidities 



Fig. ^ displays the calculated nuclear modification factor for pion production in mini- 
mum bias dAu {upper row) and in central AuAu {lower row) collisions at ^/s = 200 AGeV 
at different rapidities, t] = 0, 1, 3.2. The dAu results indicate the validity of our LO per- 
turbative QCD calculations including multiscattering and shadowing [ElE]- Data on pion 
and charge hadron production are shown from PHENIX [^j and BRAHMS 
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Figure 1. Nuclear modification factors at different rapidities: R^au '^^'^ collisions 
{upper row), RauAu in — 10% central AuAu collisions {lower row). Data are from 
PHENIX [Illlig and BRAHMS [0 d at = 200 AGeV. Full lines display our pQCD 
results, the bands indicate the theoretical uncertainties on the value of the opacity. 

In AuAu collisions strong final state interaction (induced jet energy loss) appears and 
data display suppression at high [T]. At midrapidity, our pQCD result with average 
opacity of L/X = 4.0 ± 0.25 is consistent with data at pt > 4 GeV/c, where the high-p^ 
region sets in (see Fig. H]). Recent BRAHMS data [El Ej may indicate rapidity inde- 
pendent suppression for pion production in the intermediate-pr region (p^ < 3.5 GeV/c). 
Assuming that this finding is also valid for the high-p^ region, we wish to achieve a similar 
high-pT- suppression at large rapidities. A smaller opacity, L/X = 3.5 ±0.25 at rj = 1, and 
L/X = 1.5 ± 0.25 at ?7 = 3.2 is needed to do this. (Note the difference between charge 
hadron and pion suppression.) This result indicates that longitudinally traveling partons 
see less colored matter then those traveling in the transverse direction. 

Fig. 121 displays available data and our pQCD results for AuAu collisions at ^/s = 
62.4 AGeV for different rj values. Due to the smaller Bjorken energy density [El ITij . 
e^j^j"^ < e'^j, we have used reduced opacity parameters, namely L/X = 3.25 ± 0.25 at 
mid-rapidity, L/X = 2.75 ± 0.25 at = 1. and L/X = 0.75 ± 0.25 at r/ = 3.2. Extended 
bands indicate the size of the theoretical uncertainties on the induced energy loss. 
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Figure 2. Nuclear modification factor, R\uAui at 77 = 0, 1 and 3.2 in — 10% central AuAu 
collisions at RHIC energy, ^/s = 62.4 AGeV. Our pQCD results with varying opacities 
(full lines and bands) are compared to PHENIX [ IJ, STAR [1^, and BRAHMS [i3j data. 

Comparing R^^Au results of Fig. [T] and Fig. |21 the decreasing tendency of the opacity 
parameter is clearly seen as t] is increasing. At large forward rapidities the interplay be- 
tween stronger shadowing and weaker quenching effects maintains a rapidity independent 
nuclear modification factor (no data are available at 77 = 3.2 for ^/s = 62 A AGeV AuAu). 

Whether vr suppression is approximately rapidity independent will of course be decided 
by the data. The present hint of this behavior motivated us in this study to examine 
the interplay of rapidity dependence and jet energy loss. One particle quark tomography 
combined with jet-jet correlation results provides information on the geometry of the hot 
region. Energy dependence and CuCu results may be used to verify our conclusions. 
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